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ABSTRACT 

We introduce a technique to measure gravitational lensing magnification using the variability of 
type I quasars. Quasars' variability amplitudes and luminosities are tightly correlated, on average. 
Magnification due to gravitational lensing increases the quasars' apparent luminosity, while leaving 
the variability amplitude unchanged. Therefore, the mean magnification of an ensemble of quasars 
can be measured through the mean shift in the variability-luminosity relation. As a proof of prin- 
ciple, we use this technique to measure the magnification of quasars spectroscopically identified in 
the Sloan Digital Sky Survey, due to gravitational lensing by galaxy clusters in the SDSS MaxBCG 
catalog. The Palomar-QUEST Variability Survey, reduced using the DeepSky pipeline, provides vari- 
ability data for the sources. We measure the average quasar magnification as a function of scaled 
distance {r/R2oo) from the nearest cluster; our measurements are consistent with expectations as- 
suming NEW cluster profiles, particularly after accounting for the known uncertainty in the clusters' 
centers. Variability-based lensing measurements are a valuable complement to shape-based techniques 
because their systematic errors are very different, and also because the variability measurements are 
amenable to photometric errors of a few percent and to depths seen in current wide-field surveys. 
Given the data volume expected from current and upcoming surveys, this new technique has the 
potential to be competitive with weak lensing shear measurements of large scale structure. 

Subject headings: gravitational lensing:weak ~ galaxies: active - quasars: general - galaxies: clusters - 
large-scale structure of Universe - methods:data analysis 



1. INTRODUCTION 

Gravitational lensing of an object by mass along the 
line of sight produces two measurable effects: shape dis- 
tortion and magnification of the source. It is usually im- 
possible to measure the level of magnification since one 
typically does not know the intrinsic fiux of the source. 
Because of this fundamental problem, the field of weak 
lensing has focused on shape measurements of sources 
since one can usually assume that galaxies have no pre- 
ferred orientation. (For a r eview of gravitational lensing 
theory and techniques, see [Schneider et aT]|2004 ) How- 
ever, if the intrinsic luminosity of sources can be esti- 
mated, then lensing magnification can be used to study 
large scale structure. 

Recently, large scale sky surveys have become effi- 
cient and prolific enough such that variability stud- 
ies can be conducted on large samples of rare ob- 
jects. For example, ensembles of type I (broad-line) 
quasars have been studied for the purpose of probing the 
physical causes of their common and dramatic variabil 
itv (e.g.. iVande n Berk et al.ll2004 iDe Vries et all 
IWilhite et al'l l2n08. B auer et al.ll20n9mKellv et al.ll 
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iMacLeod et"al] 120101 iMeusinger et all 1201 1[ ). In these 
works, the mean ensemble variability amplitude has been 
seen to correlate with a number of properties of the ob- 
jects: e.g., time scale of the variability, wavelength of 
observation, redshift, mass, and luminosity of the sys- 
tems. In particular, a strong anti-correlation has repeat- 
edly been observed between variability amplitude and 
quasar luminosity, particularly when comparing quasars 
of similar mass. 

Given a measurement of a quasar's variability ampli- 
tude, one can estimate its luminosity using this empirical 
relation. If the quasar were gravitationally lensed by in- 
tervening mass, the magnification effect would alter the 
observed luminosity. However, the fractional variabil- 
ity would remain unaffected, as the magnification, which 
is multiplicative on the luminosity, cancels. Therefore, 
magnification due to gravitational lensing will shift the 
quasar's position in variability-luminosity space. Quan- 
tification of this shift constitutes a measurement of the 
quasar's lensing magnification. Because the variability- 
luminosity relation has a large intrinsic variance, the es- 
timate of a single quasar's magnification is not precise. 
However, because the relation is well-determined in the 
mean, an ensemble of quasars can yield a significant mea- 
surement of the lensing effect. 

In this paper we take advantage of the variability- 
luminosity relation seen in type I quasars to measure 
their magnification due to galaxy clusters along the 
line of sight. We use the Palomar-QUEST Variabil- 
ity Survey to measure the lensing magnifications of 
3573 quasars, and compare the signal to that expected 
from gravitational lensing due to the 13,823 galaxy clus- 
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ters in the Sloan Di gital Sky Survey MaxBCG catalog 
(|Koester et all [20071 ) . We measure an average cluster 
profile shape that is consistent with previous studies of 
the MaxBCG cluster catalog, showing that this method 
of measuring lensing magnification is indeed effective. 

2. AGN STRUCTURE AND QUASAR OPTICAL 
VARIABILITY 

An active galactic nucleus (AGN) is typically described 
as having several prominent structural components. The 
central black hole is fed by an accretion disk, which 
radiates a continuum of fiux primarily in the optical 
and UV bands. Offset from the plane of the disk ex- 
ist clumps of gas which absorb the disk's flux and rera- 
diate it as broad emission lines (hence the name broad 
line region (BLR)). Emission from the disk and the BLR 
constitute the majority of the optical fiux observed in 
type I quasars, which are the highly luminous AGN for 
which the observer has a direct line of sight to these two 
structures. (Type I quasars will hereafter be referred to 
simply as quasars.) The detailed geometry of the disk 
and BLR are not well unders tood; they m ay be p art of 
one continuous structure (e. g. iNicastro et al. 20031, lElvisI 
120001 [Lovegrove et al.l[2010( ). but they traditionally are 
described as being distinct. A minority of quasars also 
show evidence of radio emission from a jet. The fraction 
of quasars that have strong radio emission (i.e. are ra- 
dio loud ) is uncertain; esti mates range from roughly 5 
to 25% ( Zamfir et al7'2008h. Jet emission can account 
for a portion of a quasar's optical flux, but the fraction 
of optical flux that is jet-based depends strongly on the 
properties and orientation of the particular object. 

Radio-loud quasars are on average more optically vari- 
able than radio-quiet quasars; in particular, their short- 
timescale variability is enhanced (e.g., Kell y et al. 200!|). 
This distinction implies that the jet contributes to the 
optical variability in these systems, most likely due to 
the dissipation of shocks. 

Not all quasar variability is due to jet physics; quasars 
with no ev idence of a j et als o exhibit significant flux 
variations. ISesar et al.l (pOOl . using the Sloan Digital 
Sky Survey (SDSS), found that at least 90% of quasars 
are variable by at least 0.03 magnitudes on timescales 
of up to several years. The variability of individual 
nearby radio-quiet quasars (and their lower-luminosity 
counterparts, Seyfert I galaxies) have been studied in 
detail using the technique of reverberation m apping (e.g. 
iBlandford fc McKedlTOSl iBentz et al.ll2009f ): spectra of 
the AGN are regularly obtained, and the variability of 
the spectral features is measured as a function of time. 
Such monitoring has revealed that the disk continuum 
fluctuations are followed by fluctuations in the broad 
emission lines that echo the continuum variation. The 
time lag between the continuum and broad line variabil- 
ity constitutes a measurement of the BLR size, which is 
typically on the order of light-days. These measurements 
show that the primary source of variability in these ob- 
jects is inside the BLR, presumably from the accretion 
disk and instabilities therein. 

Accretion disk instabilities have long been a prime 
suspect as the main source of quasar variability. 
iKawaguchi et al.l ()1998[ ) calculated the expected rela- 
tionship between variability amplitude and timescale 
based on a model of disk instabilities, as well as a 



relation expected fr om s tarbursts in the quasar host 
galaxy. iHawkinsI (|2002[ ) extended this work to in- 
clude a prediction based on microlensing of the quasars. 
A number of studies of ensemble quasar variability 
have measured the average variability amplitude ver- 
sus timescale relation, which best agrees with the ac- 
cretion dM^JristaM & Pctcrsop 
120011 iVanden Berk et all 12004 [Pe Vrics ct aL ■200E 



Wilhite et al.l 120081 iReiigstorf et alT 2006. B auer et alj 
l2009bl iMeusinger et al.l l2011[ ). However, this type of 
analysis would benefit from updated predictions using 
more complex models. 

Disk instability models invoke relatively small, short- 
lived flares from numerous instabilities that combine to 
produce a stochastic light curve. These fluctuations 
could be prompted by changes in, for example, the ac- 
cretio n rate (e.g.. | Li fc C ao 2008) or in a magnetic fleld 
( Hiroseend1l2009( ). Stochastic light curves have been 
succesfully used to d escribe optical observations of large 
samp l es of quasars (iKellv et all 120091 iKoztowski et al.l 
120101 iMacLeod etal . 2010). These works model the 
quasar light curves as a damped random walk, using only 
three parameters: the typical amplitude of the short- 
timescale variability, the damping timescale, and the av- 
erage flux of the light curve. Such light curves, generated 
from a constant component plus similar flares with ran- 
dom start time, follow the variability-luminosity relation 
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where L is the object's l uminosity and 5 depends 
on the details of th e flare s (jCid Fernandes et al.l Il996l 
iCid Fernandes et"al1 120001 ). Qualitatively, the anti- 
correlation simply notes that, given similar flare prop- 
erties between two quasars, the fractional variability 
/S.L/L will be smaller in the quasar whose mean lumi- 
nosity L is larger. This prediction indeed describes the 
variability-luminosi ty anti-co rrelatio n seen in quasar e n- 
sembles (Vandcn B erk et al. l 2004, Ba uer et al.ll2009b[ ). 

Alternatively, it is possible that much of the observed 
relationship between variability and luminosity is a sec- 
ondary effect of a correlation between, for example, vari- 
ability and the quasar's Eddington ratio. The Edding- 
ton ratio of a quasar describes its accretion rate, and 
can be estimated as a ratio be tween the quasar's lumi- 
nosity and its mass. lAi et al.l (2010) flnd the correla- 
tion between Eddington ratio and variability to be more 
significant than that between luminosity and variability, 
using a sample of several hundred SDSS Seyfert I galax- 
ies. The authors note the positive correlation between 
the accretion rate and the disk radius which emits at a 
given wavelength. Therefore an anti-correlation between 
variability amplitude and Eddington ratio implies that 
the outer regions of the disk show less variability than 
the inner radii (regardless of the physical mechanism be- 
hind the variability). This position-dependence qualita- 
tively agrees with the observed anti-correlation between 
variability and wavelength (Wilhite et al.l 120051 ) , if the 
hotter, bluer, inner areas of the disk vary more strongly 
than the redder, outer regions. The variability amplitude 
may be anti-correlated with the quasar's luminosity via 
such a positional dependence, and therefore through the 
definition of the Eddington ratio as related to the lumi- 
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nosity. As described in section[31 we normalize the quasar 
variability measurements according to the objects' phys- 
ical properties, including black hole mass. This proce- 
dure effectively allows us to measure how the variability 
amplitude scales with luminosity, for quasars of similar 
mass. Holding the mass constant in this way allows the 
variability-luminosity relation to contain the same infor- 
mation as the correlation between the variability and the 
Eddington ratio. 
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3. MEASURING ^l 

The average amplitude of quasar variability has 
been seen to depend on several factors: time lag 
between measurem en ts, l uminosity of the quasar 
(iVanden Berk et al l [200l iDc Vric s et al.l l2005l 
Wilhite et al.l 120081 iBauer et al.l i2009b. iMacLeod et al.l 
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2011) 



iMeusineer et ah. 

serva ti on ("Vanden Berk et al 
[2005L iMe usinecr et al.l I2011D 



wav el ength of ob- 
I200I iDe Vries et all 

^ . and bla ck hole mass 

(jWilhite et al. " 20081 lBauer"etaII l2009bD . The depen- 
dence of variability a mplit ude on redshift is less obvious; 
iVanden Berk et ahl (|2004D measured a slight increase 
in the variability with redshift, while iDe Vries et al.l 
(I2005D measure d a sl ight de crease. More re cently , 
IMacLeod et"all (|2010f ) and iMeusingeteLaU (|201lD 
have measured no significant dependence of variability 
amplitude on redshift. 

In practice, the variability-luminosity trend measured 
is of the form: 



log{V) = C -ax log(L,„oas)- 



(2) 



A linear relation has indeed be e n ob served 
(IVanden Berk et al.l I200I IBauer et al.l l2009bD . al- 
though there are conflicting results for the value of the 
power-law slope a, perhaps due to selection effects. If 
faint quasars are included in the analysis for which one 
cannot observe the full extent of the variability, the 
measured slope will become artificially shallow; this 
effect most clearly manifests itself as a flattening in 
the relation at the lowest observable l uminosities, as 
is illustrated in figure 5 of IBauer et all C2009b,) . The 
low-luminosity limit of the binning scheme in this work 
is chosen to exclude this regime from the data set. The 
value of the constant C depends on the details of the 
normalization of the data, as described below. 

When studying how the variability of a large quasar 
sample depends on one of the quasars' properties, we 
must treat the parameters as independently as possi- 
ble. To this end , we use a method in troduced by 
IVanden Berk et all (2004^ and adopted in I Bauer et all 
(|2009b( ). Four basic quantities are known for all of the 
quasars in our sample: time lag between measurements 
r, luminosity L, estimated black hole mass M, and red- 
shift z. There are known correlations between all of these 
parameters, due to physical relationships or artificial ef- 
fects such as detection biases in flux-limited surveys. To 
avoid these complications and study only the dependence 
of variability on luminosity, we would like to identify a 
set of quasars with identical properties except for their 
luminosity, and then examine how the variability differs 
between them. To approximate this procedure, we have 
split each parameter's range into bins: 8 bins in r, 6 
bins in Af, 6 bins in z, and 4 bins in L. The bin limits 



Table 1 

Bin limits used in determining the normalization constants. 
Measurements with values outside the limits are not considered. 
Units of time lag t: days; black hole mass M: 10* X Mq, 
luminosity L: ^'j^^ . 



are given in table |TJ quasars with properties outside the 
given ranges are not used in the analysis. 

To measure quasar variability we use a quantity similar 
to that of the structure function. We define: 



V = v/(Am)2-a2, 



(3) 



where Am is the magnitude difference between two in- 
dependent observations of an object, and a is the error 
on those measurements. This is similar to the s truc- 
ture functio n as use d in IVanden Berk et al.l ()2004D and 
IBauer et al.l (|2009bD : however, here instead of being an 
ensemble measurement, one V is measured for each pair 
of magnitude measurements of a quasar. Four mea- 
surements of a quasar will yield six Am measurements, 
and therefore six different V measurements for the sin- 
gle quasar. As V is imaginary when Am is less than 
the measurement error cr, we only use data which show 
significant (> la) variability. This cut on the data is 
described further in section [5.11 

For each multi-dimensional bin, a mean variability am- 
plitude V is determined by taking the mean of all V val- 
ues measured for that bin. Then, holding constant the 
indices for time lag, mass, redshift, and wavelength, one 
can compare the V values across the 4 L bins. This pro- 
cedure yields 8 x 6 x 6 = 288 possible 4-point plots of 
mean variability amplitude V versus luminosity, or 1152 
possible V values. Most of the multi-dimensional bins 
are not well populated by the quasar sample (for exam- 
ple, high- redshift low-luminosity bins). In fact, we obtain 
403 bins with at least 50 measurement pairs, which is the 
minimum we require in order to adequately determine V. 
To examine the overall behavior of V with respect to L 
one can normalize the 4-point measured trends together 
and average the resulting normalized data in each L bin 
to find a simple, meaningful result of how the variabil- 
ity scales with the quasar luminosity. The normalization 
consists of an additive constant in log(V^), i.e. each 4- 
point measured trend has its own constant C as defined 
in equation |2j The t, M, and z multi-dimensional bin 
that has the best statistics is chosen to be the standard, 
and the log(T^) versus log(L) trends from all other r, 
M, and z bins are normalized to that standard using 
one constant offset per r, M, z combination. The con- 
stant is determined by minimizing the chi square differ- 
ence between the V values from the two datasets in the 
same L bin, for the L bins where there exist data from 
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both sets. For a visual repres entation of the norm aliza- 
tion procedure, see figure 4 in IBauer et alJ (I2009bn . Af- 
ter averaging the normalized data, we are left with one 
Kiorm versus L trend with arbitrary y axis normaliza- 
tion but meaningful slope. This normalization technique 
has been shown to give results for variability amplitude 
versus time lag r that are co nsistent with i ndepen dent 
measurements (see table 4 in IBauer et al.l ()2009b[) . In 
this way, we study how the variability scales with lumi- 
nosity, comparing only objects with similar values of the 
other parameters. After the normalization, deviations 
from the mean Viiorm — L relation will not be caused 
by known, but lensing-independent, correlations such as 
that between variability and time lag. Using the normal- 
ization constants calculated in this way, each measured 
variability amplitude V is normalized according to its r, 
M, and z bin; the resulting V,iorm then can be used to 
estimate the quasar's lensing magnification, as detailed 
below. 

We note that, for data taken in a single pass-band, 
the redshift of the quasar is degenerate with the rest- 
frame wavelength observed: the effective wavelength of 
the filter in the quasar's rest frame scales as (1 -I- z). Be- 
cause rest-frame wavelength of observation scales mono- 
tonically with redshift, these two variables are simultane- 
ously accounted for in our normalization versus redshift. 

When fitting for a and C in equation [2] , we find the 
linear fit that best describes the whole quasar sample. 
On average, the measured luminosity Lmoas of quasars in 
the sample is equal to /I-^intiinsic, where /I is a magnifi- 
cation value representative of the entire sample. We can 
therefore rewrite equation [5] as 



log(Korm) = C - a X l0g(/lLi, 



c) 



(4) 



where we have also explicitly noted the fact that the V 
measurements are normalized. 

Once we have determined the unique constants a and 
C for the whole data set, we can use a measured and 
normalized l^orm of a quasar to constrain /J^intrinsic since 
equation m is equivalent to 



10^ 



l/a 



(5) 



Qualitatively, we simply use the Vnorm for a quasar to 
read off an expected value of Lmeas from the linear 

l0g(l^„orm) - log(L„icas) plot. This CXpCCtcd Lmcas, aS 

it is the prediction calculated using the luminosities of 
the entire quasar sample, is an estimate of /I-^intrinsic, 
where ^intrinsic is the intrinsic luminosity of that quasar. 
Comparing this value to a direct measurement of the lu- 
minosity of the quasar Lmoas, we can measure any shift 
in the luminosity, i.e. magnification or demagnification: 



Mtrue 



I/O 



(6) 



We define /i' = /Xtruo//Jj the relative magnification with 
respect to the typical value for the whole sample. For a 
large set of quasars taken from a broad sky area, such 
as the data used in this work, we expect JI to be close 
to unity. In this way, each variability measurement V, 
calculated from a pair of magnitude measurements Am, 
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Figure 1. RG-610 filter transmission, multiplied by CCD re- 
sponse. 

yields an estimate of the quasar's relative magnification 

4. PALOMAR-QUEST RG-610 DATA SET 

We measure the lensing magnification of quasars us- 
ing the Palomar-QUEST Variability Surv ey. The sur- 
vey u sed the QUESTII large area camera ()Baltav et al.l 
I2OOI . which comprises 112 CCDs with a pixel scale of 
0.878" /pixel, with an overall field of view of 9.6 square 
degrees. The camera was mounted on the 48" Samuel 
Oschin Schmidt telescope at the Palomar Observatory, 
during which it aquired two main data sets. First com- 
pleted was a 3.5 year multi-color survey covering 15,000 
square degrees roughly four times in each of seven opti- 
cal filters (Johnson UBRI and SDSS riz), with a depth 
of roughly mag 21 in the SDSS r filter. The se data have 
been used to study the var iability of AGNs () Bauer et al.l 
'2009bl lBauer et al.ll2009al) . to searc h for highly variable 
objects in both the archival dat a (IBauer et al] l2009cD 
and real-time d ata (jP iorgovski ct al. 2008' ), and to study 
brown dwarfs (jSlesnicb i2006b . .Slesnick i2006aD . among 
other phenomena. The multi-color survey was executed 
concurrently with a 5 year, 30,000 square degree sur- 
vey which used a single broad, red, RG-61G filter, and 
reached a depth of about mag 19.5 in SDSS r in each 
60 second exposure (under dark sky conditions). The 
bandpass for the RG-610 filter (filter transmission times 
the QUEST CCD quantum efficiency) is shown in figure 
[TJ The single-filter survey has been used a s the discov- 
ery d ata for the Nearby Supernova Factory (jCopin et al.l 
|2006[ ). as well as for searches for dwarf p lanets and other 
solar system bodies ()Brown et al.ll20Q5f l. 

This work uses the single-filter Palomar-QUEST data, 
and constitutes the first results from these data which 
rely on precision photometry. For this reason, we de- 
scribe the DeepSky image reduction pipeline and the 
photometric calibration techniques which are used to 
generate flux measurements. 

4.1. The DeepSky Image Reduction Pipeline 

The Palomar-QUEST RG-610 data were taken by a 
variety of groups with different strategies, under differ- 
ent observing conditions, and with no uniform calibra- 
tion procedure; daily operations in each of the groups 
emphasized simply discovering transients, rather than 
more precise photometry. Industry-standard reduction 
frames such as dark frames, dome flats, twilight fiats, and 
fringe frames, were taken rarely or not at all. The Deep- 
Sky project has re-reduced these data using a standard 
pipeline designed to improve the precision of differential 
photometry, which we briefly describe below. 
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The DeepSky reduction suite treats each of the 112 
CCDs in the QUESTII camera as a separate detector, 
attempting only to produce a uniform response on each 
frame rather than trying to flatten the detector response 
across the entire mosaic. Data frames for each stage 
of the reduction (dark, flat, etc.) were pooled for each 
month, and sometimes from several months in succes- 
sion depending on what was available, to achieve the high 
statistics necessary to reject bad pixels and regions of the 
camera with unstable performance. The algorithms were 
cross-validated on subsamples of the data for months 
with many reduction frames, suggesting that the loss of 
performance is minimal when building reduction frames 
monthly rather than more frequently (weekly or daily). 
Each month about 100 calibration images from each chip 
were used to create 'superframes' for the reduction. 

Bias levels were first subtracted from each 600 x 2400 
pixel image using a 40 x 2400 overscan region. Dark 
frames of different exposure times (10, 60, 100, and 
240 sec) were combined, with the dark current Di at 
each pixel i being fit to a model linear in the exposure 
time t: Di = ai + bit, where the 'instantaneous' portion 
corresponds to dark current which accumulates during 
readout. Iterations of the fitting procedure rejected out- 
lier values in determining the model parameters ai and 
bi. Pixels for which no solution converged, or for which 
the "1 cr" interquantile width (median minus 16th per- 
centile) of the distribution of residuals from the fit was 
more than three times wider than the median residual 
value for the image, were masked as bad pixels. All im- 
ages were then corrected according to their exposure time 
using the superdark Di. 

In general neither twilight nor dome flats were regu- 
larly available, and images taken under dark-sky condi- 
tions were often contaminated with strong fringes. We 
therefore chose to construct flat fields using science im- 
ages taken bot h during twilight and i n moo nlight. The 
calculations of iKrisciunas fc Schaefeii (|1991[ ). appropri- 
ately adapted to Mt. Palomar, suggests that the loss of 
flatness due to the non-uniformity of moonlight should 
be less than 2% over the extent of each chip (0.5 de- 
grees) for images taken more than 20 degrees from the 
moon, and we required this of images used to build the 
flat field. We also required that the sky level from scat- 
tered moonlight in these images be at least five times 
that typical of dark sky conditions. Images passing both 
of these cuts were normalized to a mean sky brightness 
of 1.0, and combined by taking the median at each pixel 
location to produce a superflat for the month. Pixels 
with unusually large fluctuations relative to other pixels 
were masked as bad, as for superdarks above. All images 
were then flattened using these superflats. 

Superfringes were constructed in an analogous man- 
ner to the superflats using images taken under dark-sky 
conditions. The baseline fringe pattern was assumed 
to be a time-invariant characteristic of each chip over 
the month. The fringe amplitude in each science image 
was determined through cross-correlation with the super- 
fringe, and the fringes in the image removed by subtract- 
ing the superfringe scaled by this amplitude. Changing 
sky conditions can of course produce time variations in 
the fringes, particularly in such a wide-band filter. We 
accept this as a systematic error in the photometry, with 
a contribution roughly on the same order as the sky noise 



under dark sky conditions. 

Object detection and aperture photometry were 
performed on the detrended (bias-subtracted, dark- 
subtracted, flattened and de-fring ed) images using 
SExtractor (Berli n fc Arnoutslll996( ). Astrometry was 
perfor med on each fr ame using the astrometry.net 
suite (jLang et al.ll2009[ ). Magnitudes were measured for 
each object using ten aperture diameters ranging from 
1 pixel to 16 pixels; this work uses the 3-pixel diam- 
eter aperture measurements as the primary flux mea- 
surement. Aperture corrections were calculated as the 
clipped median of the 12-pixel diameter aperture flux 
divided by the 3-pixel aperture flux. This correction 
was calculated using all good-quality measurements on 
a frame, and calculated separately for each frame. 

4.2. Photometric Calibration 

A number of cuts are made on the quality of the 
data. First, SExtractor flags are checked for objects 
that are blended or saturated, close to image bound- 
aries, or for which the measurement failed; these flags 
eliminate about 3% of the data. A further saturation 
check is made by comparing object fluxes to saturation 
levels determined for each chip, which are more accu- 
rate than the single baseline saturation level checked by 
SExtractor. Very few objects are eliminated by this 
cut. Variations in the background sky level on the small 
spatial scale of 5x5 pixels are examined, and a frame is 
discarded if they are typically large enough to contribute 
a 3% error to the photometry of an object with magni- 
tude 18; this removes about 3% of the images. Frames 
are also rejected for which the median full width at half 
maximum (FWHM) exceeds 5 pixels, eliminating about 
6% of the data. To ensure that closely spaced detec- 
tions do not contaminate the aperture photometry, an 
object is rejected if it has a neighbor within four arc sec- 
onds (~3.5 pixels); most object affected by this cut are 
already discarded for being blended in SExtractor. 

Spatial fluctuations of the moon/twilight flat fields and 
variability in night sky lines during observations made it 
difficult to completely remove the fringes from the data. 
This results in a ~ 2% fiux systematic for bright objects, 
which is added in quadrature to the statistical errors. 

A standard method of calibration is simply to de- 
termine a multiplicative constant zero point which cor- 
rects differences in average response (both intrinsic and 
weather-related) between two images. We perform such 
a frame-based calibration as a first step. All overlapping 
images are determined, and the one in which objects have 
the brightest measurements is assigned as the standard. 
A single zero point is determined for each frame in or- 
der to bring it to the level of the standard. This con- 
stant is determined using objects of all magnitudes. The 
correction is also calculated separately for bins of differ- 
ent magnitudes: five bins between SDSS r magnitudes 
of roughly 15 and 20. On average, the difference in the 
mean correction between the brightest bin and a fainter 
one ranges from about 3% to 6%. If a larger variation 
than this is seen, the frame is assumed to be non-linear 
and is not used. This eliminates about 15% of the data; 
the non-linearity is partially due to intrinsic properties 
of the CCDs, but primarily due to the insufficient quality 
of the fringe and flat fleld calibration data. Poor frames 
are also identified by comparing their measurements to 
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the mean measurements of the objects they hold. If 15% 
of the flux measurements on a frame disagree with their 
objects' means by more than 3 times the measurement 
errors, then the frame is discarded. This removes about 
8% of the data, for which the systematic effects are catas- 
trophic. 

The systematic residuals in the data make a frame- 
based calibration insufficient for accurate photometry; 
while a systematic term describes the typical uncertain- 
ties, the error distribution is not Gaussian and there- 
fore this procedure generates many photometric outliers. 
Because the spatial scale of the systematic variations is 
much larger than the size of each object, we can refine 
the frame-based calibration using an object-based one. 
This technique involves generating a zero point for one 
object at a time by considering only its close neighbors. 
A large scale systematic variation far away from an ob- 
ject will affect a frame-based calibration; an object-based 
one will be less sensitive to such features. We therefore 
implement, as a second step, an object-based calibration. 
We use the average of at least 10 neighbors within 150 arc 
seconds to calculate the correction between an object's 
measurements on two scans. About 8% of the detections 
do not have enough neighbors to satisfy this requirement, 
and are therefore disregarded. 

The frame calibration is calculated using only ob- 
jects detected by the Sloan Digital Sky Survey (SDSS; 
lAbazajian et al.H2009f ) and classified by them as point- 
like. This precaution eliminates the possibility of our 
inadvertently using data artifacts, or galaxies which may 
exceed the aperture size, as calibration objects. For 
the purposes of this work, the object calibration is only 
performed on quasars which have been spectroscopically 
identified by SDSS. The reference objects used in the ob- 
ject calibration are limited to pointlike SDSS detections. 

The strict cuts implemented in the photometric cali- 
bration throw away roughly 35% of the measurements. 
The systematic error is typically 2%, with about 3% of 
measurements disagreeing with the mean flux of the ob- 
ject by over 3 times their measurement errors. From 
previous variability studie s we expect < ~1% of m easure- 
ment s to be truly variable (jHuber et al.l[2006; .Sesar et al.l 
12001 : therefore, we infer that our photometry contains 
residual calibration errors in ~2% of the measurements. 
Because in this work we deal only with the ensemble av- 
erage over hundreds of measurements, and we expect the 
systematic error to be uncorrelated across different ob- 
servations, we do not expect the residual errors to affect 
our results. Nonetheless, as a precaution against these 
errors we make additional cuts against variability behav- 
ior uncharacteristic of quasars; this is described further 
in section [5?T] 

5. RESULTS 

We study a sample of 3573 spectroscopically identified 
quasars, spanning roughly 5,000 square degrees, from the 
SDSS data release 7. For each quasar we calculate its 
luminosity at 2500 A by multiplying a redshifte d com- 
posite quasar spectrum taken from Vandcn Ber k et al.l 
(|200l with SDSS fiher curves. We use the SDSS z 
band magnitud e, corrected for galactic extinction using 
the maps from iSchlegel et al.l |l998) , to normalize the 
composite spectrum's amplitude. We choose the z band 
as it is least extincted by dust along the line of sight. 



We then take the normalized fiux at 2500 A and convert 
it to luminosity using the object's spectroscopic redshift 
and the cosmological parameters, which we assume to be 
1] = = 0.7, = 0.3, iJo = 0.71^. We also 

estimate the black hole mass for each quasar from the 
M — a relation, usi ng the widths of the H /^ and Mgn 
emission lines as in iSalviander et al.l ()2007[ ) , and calcu- 
lating the radius of the BLR as described in iBentz et al.l 
(i2009). 

In many cases during the Palomar-QUEST Survey, the 
same sky area was observed twice in the same night. Be- 
cause quasars are not variable on such short timescales 
to the precision of our measurements, we average all in- 
tranight measurements to yield one flux value per night. 

5.1. Data Cuts 

As described above, we have made restrictions on the 
data quality in order to produce reliable photometry. 
Here we describe the criteria we use to select our quasar 
sample so that we have a well-measured, homogeneous 
dataset. 

While the majority of quasar variability is thought to 
be due to accretion disk instabilities, some AGN show 
blazar characteristics, with considerable flux and vari- 
ability from a jet component. Jet-based variability is 
much more dramatic than that seen in typical quasars, 
and because it is powered by shock dissipation and en- 
hanced by relativistic beaming, it is not expected to 
follow the variability-luminosity relation seen in most 
quasars. It is therefore important to restrict the quasar 
list to those displaying variability that is characteristic 
of quasars. 

We select quasars from SDSS that have the following 
properties: 

1. spectroscopic identification in the SDSS, with a 
match to the spectral cross-correlation template 
number 29 (QSO) with confidence of at least 0.95 

2. no det ection in the FIRST JBecker et~allll995[ ) or 
RASS ([Anderson et al.ll2007[ ) catalogs 

3. redshift between 0.4 and 2.2 

4. estimated black hole mass between lO^M© and 
7.5 X lO^Mo 

5. luminosity between lO^^-^'^fjl and 10^^'^ j^. 

Criteria 1 and 2 are meant to select radio-quiet AGN 
whose optical fiuxes do not show significant contribu- 
tions from a jet. Numbers 3-5 serve in part to crop the 
tails of the parameter distributions, and are the edges 
of the binning scheme that we use when normalizing 
the data as described below. The bin limits are stated 
again in table [T] The low redshift cut is chosen to elim- 
inate objects that may appear spatially extended in the 
data, and therefore be poorly measured by the Palomar- 
QUEST photometry. The high redshift cut ensures that 
the quasar sample is selected homogeneously using UV- 
excess techniques. The low luminosity cut is chosen be- 
cause for objects fainter than this level we see a flattening 
in the variability-luminosity relation, implying that these 
less luminous quasars vary to fluxes below our detection 
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limit, causing us to not observe their full range of vari- 
ability. 

These criteria yield a sample of 4845 quasars that have 
Palomar-QUEST repeated measurements. 

Furthermore, the quasars' Palomar-QUEST measure- 
ments must satisfy: 

1. V < 0.1 where V is defined as in equation [3] and is 
calculated using all measurement pairs with rest- 
frame time lag less than 10 days 

2. V < 0.5 where V is defined as in equation [3] and is 
calculated using all measurement pairs with rest- 
frame time lag greater than 100 days. 

These cut values are motivated by the va riability am- 
plitude s of quasars and blazars measured in I Bauer et al.l 
(|2009af ). Because the variability amplitude is unchanged 
for lensed and unlensed objects, this is an unbiased cut 
with respect to lensing analyses. This step eliminates 967 
objects, leaving 3878 quasars. We note that this cut elim- 
inates a much larger fraction of the AGN than is expected 
from the relative numbers of known blazars and quasars. 
However, the cut is not only sensitive to jet-based vari- 
ability but also calibration errors, and is a strict but im- 
portant criterion to ensure that poor quality measure- 
ments do not overwhelm the lensing signal. These 3878 
quasars have, in total, 230,674 Palomar-QUEST mea- 
surement pairs. 

We only use measurement pairs for which: 

1. the rest frame time lag between measurements is 
shorter than 400 days 

2. the measurements show significant (> la) variabil- 
ity 

3. the variability amplitude is less than 1 magnitude 

4. no more than 40 measurement pairs from each 
quasar are included. 

Cut 1 eliminates a long but low-level tail out to larger 
time lags, and throws away 19,523 pairs. Cut 2 is nec- 
essary to avoid the calculation of imaginary V values 
for individual measurement pairs and is very restrictive 
because the Palomar-QUEST systematic error is of the 
order of typical quasar variability amplitudes over rest 
frame timescales of weeks. Still, as will be seen in figure 
m the data after this cut continue to follow a variability- 
luminosity relation described by equation [51 and there- 
fore can be used to measure magnification using equation 
[HI We also note that removing the lowest signal-to-noise 
V values reduces error due to imprecise determination of 
the survey's systematic uncertainties, given the depen- 
dence of V on the measurement uncertainty as stated in 
equation [31 Because the systematic uncertainties stem 
in part from imperfect flat field and fringe calibration 
images, they can be time and position dependent and 
therefore difficult to determine precisely for each mea- 
surement. This cut eliminates 134,745 pairs. 

Cut 3 is a final check against outliers that are unchar- 
acteristic of quasar behavior, and only throws away 17 
measurement pairs. Cut 4 serves to keep a minority of ex- 
ceptionally well-measured quasars from dominating the 




log( L2500 ) 

Figure 2. Log Variability vs. Log Luminosity at 2500A, after 
normalizations, determined using 57,359 measurement pairs. The 
best-fit line is shown, with slope —0.565 it 0.007. 

dataset, since the number of measurement pairs rises as 
the number of measurements squared. This cut removes 
14,671 pairs. 

Normalization of the quasars, described in section [3l 
is only done if there are sufficient quasars with similar 
properties so that we can determine accurate zero points. 
We require 50 measurement pairs per t M z L bin, which 
removes 4359 measurement pairs. 

We are left with 57,359 useable measurement pairs, 
from 3573 quasars. We note that many of our cuts are 
due to properties of the data such as statistics and mea- 
surement errors, rather than properties of the quasars. 
Cuts on future quasar data sets may not need to be so 
restrictive. 

5.2. Variability- Luminosity Relation 

After normalizing the data as described in section [31 
we plot the normalized variability vs. luminosity rela- 
tion, which is shown in figure [2l The y axis error bars 
on the points are errors on the mean calculated from the 
spread of the measurements in each bin. The bins in 
the figure are smaller than those used in the normaliza- 
tion, in order to show the relation in more detail. The 
trend is indeed well-described by a power law, with in- 
dex a = 0.565 ± 0.007. The best fit line is shown on the 
plot. This slope is steeper than that found in other works 
(e.g.. IVanden Berk et al .''2004. B auer et al.ll2009bf ) : how- 
ever, it is not directly comparable since we only consider 
measurements for which we observe significant variabil- 
ity- 

The precision of our magnification measurements relies 
on a narrow spread in the variability-luminosity relation. 
The residuals in log(y) around the best-fit line are dis- 
tributed normally, with a standard deviation a = 0.22 
that is constant across our range of luminosity. These 
residuals are larger than the typical measurement errors, 
and therefore reflect an inherent spread in the quasar 
variability properties. There are several possible sources 



8 



of the enhanced variance. 

The measurement error on the luminosity is small, and 
is based simply on the uncertainties in the SDSS broad- 
band magnitudes. This error estimate ignores the fact 
that the quasars vary by more than these uncertainties; 
the true mean luminosity could be tens of percent differ- 
ent from what was measured when the SDSS made its 
photometric observation. Since the spectroscopic quasar 
sample is much brighter than the SDSS photometric de- 
tection limit, there should be no Malmquist bias in the 
quasar sample and the error on the measured luminosity 
should be unbiased. Therefore, when discussing the vari- 
ance of the relation, we assume that the luminosity is cor- 
rect and simply absorb its error into the y axis spread in 
values. In future studies this uncertainty can be reduced 
in cases where there are more photometric measurements 
per quasar, which can be averaged to produce a more 
accurate typical luminosity. (This is not possible with 
the Palomar-QUEST data, as it is relatively calibrated 
rather than absolutely.) Or, if many measurements are 
available, a lightcurve may be fit to each quasar's data 
to determine the object's baseline luminosity most accu- 
rately (see, e.g., [MacLeod ct al. 2010). 

If the spread in the variability-luminosity relation is 
due to an effect with zero bias, then all quasars will in- 
trinsically lie on the mean relation and increased statis- 
tics for a given object will move it closer to the mean 
trend. If quasars in the sample have heterogeneous vari- 
ability properties, an object may lie intrinsically off the 
variability-luminosity relation and more statistics will 
only reveal this property more clearly. Some insight 
into this issue is given in iMacLeod et al.l (|2008D . which 
explores the implications of using many, as opposed to 
simply two, measurements of each quasar when studying 
the objects' variability properties. They note that using 
only two measurements yields ensemble variability re- 
sults similar to those determined using many epochs from 
the SDSS stripe 82 region. However, they also point out 
that there is an intrinsic spread in the quasars' variability 
behavior that is not currently understood. This intrinsic 
difference between individual objects surely contributes 
to the scatter observed in the variability-luminosity re- 
lation, and is our motivation for including no more than 
40 measurement pairs from each object in this analysis. 
More insight into the physical processes driving the dif- 
ferent variability properties of individual quasars would 
help in understanding and reducing the systematic errors 
in the lensing magnification measurement. 

Assuming statistic ally identical quasar li ghtcurves with 
power law spectra, iBauer et al.l ()2009b[ ) used Monte 
Carlo simulations to show that using few measurements 
per quasar to calculate V versus r gives a slope that is un- 
biased, although windowing effects can make the trend 
non-linear at very short and long time lags (which are 
excluded from the data set in this work). Therefore we 
expect no significant bias to be introduced in this anal- 
ysis due to the behavior of V in the presence of sparse 
data sampling. 

5.3. Magnification Measurements 

Now that we have calibrated the quasar measurements 
and confirmed the variability-luminosity relation, we cal- 
culate the magnification /i' = ^t-mc/'P as defined in equa- 
tion |6l The distribution of log(/i'), with one entry per 
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Figure 3. Logarithm of the measured quasar magnifications, de- 
termined using 57,359 measurement pairs. 

quasar measurement pair, is shown in figure [3l The dis- 
tribution is wide and close to Gaussian. This is due to the 
distribution of the residuals of the V — L relation, which 
has a similar shape. Because the scatter in the V — L 
relation is large, the width of the \og{^') distribution is 
much wider than the range of magnifications expected in 
the data. Just as the mean of many variability measure- 
ments produces a precise variability-luminosity relation, 
the mean of many magnification measurements produces 
a meaningful result, with an error on the mean much 
smaller than the width of the distribution. 

In order to obtain mean magnifications, we must bin 
the measurements in a physical quantity which we be- 
lieve to be related to the magnification. Here we choose 
to bin the measurements in terms of their distance per- 
pendicular to the line of sight (scaled as distance/i?2oo) 
from the nearest member of the MaxBCG galaxy cluster 
catalog. 

5.4. The Expected Signal 

The strongest magnification of quasars is expected to 
be due to clusters of galaxies along the line of sight. 
To estimate the level of magnification that we expect 
quasars to experience, we can calculate the magnifica- 
tion due to the members of known cluster catalogs. For 
this purpose we use the MaxBCG catalog, which con- 
tains 13,823 clusters and covers about 7,500 square de- 
grees. The MaxBCG catalog is estimated to be 90% pure 
and 85% complete for clusters between re dshifts 0.1 and 
0.3 an d with masses greater than IO^^Mq (jKoester et al.l 
|2007[ ). By adopting a model for each cluster, we can cal- 
culate the expected magnification of each quasar due to 
the galaxy clusters in the catalog. 

We assign a Na varro-Frenk- White (NEW) profile 
(|Navarro et al.lllQQTh to each cluster: 

^^''^ = (r/r,)[l + (r/r,)]2 
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where c is the concentration parameter of the profile, 

Sc = ^ in(i+c)-c/(i+c) ' = R200/C, and i?200 is the 
radius inside which the density is 200 xpcrit, the critical 
density of the universe at the redshift of the cluster. The 
total mass inside i?200 is -A^an n- We assume the cluster 
mass- richness relation used by IRozo et al.l ()2009( ): 

to estimate M200 given the richness A^200 stated in the 
cluster catalog. M200 and c have been seen to correlate 
with each other; we use the mass-concentration relation 
from IDuffv et al.l ([2008.1 : 

/ M \ -0.084 

c(M»„,.) = 5.71x(5-^) 

to determine c given M200 and the cluster redshift z. 
When calculating the magnifications, we truncate each 
profile at 3 X i?200- If a quasar is close to more than 
one cluster, we multiply the magnifications due to each 
nearby cluster to produce a magnification estimate that 
includes contributions from all known clusters. Multi- 
ple magnifications are a minor effect; only 6.5% of the 
quasars have expected magnifications of > 1% from more 
than one cluster. 

A histogram of the logarithm of the predicted magni- 
fications for the quasar sample is shown in figure ID We 
predict one magnification per quasar, but we measure one 
fi' per pair of magnitude measurements. The thick-lined 
histogram in figure S] shows the distribution of magnifica- 
tion predictions, with one entry per quasar and therefore 
3,573 entries. The thin-lined histogram includes one en- 
try for each magnitude measurement pair: if a quasar 
has N measurements, then that quasar's predicted mag- 
nification will have N{N — l)/2 entries of identical value. 
The thin histogram therefore has 57,359 entries and is 
comparable to figure [31 The vast majority of the quasars 
are not significantly magnified by the clusters; however, 
there is a small tail of highly magnified objects up to 

~ 2. The striking difference in the shapes of the ex- 
pected magnification distribution (figure 2]) and the mea- 
sured distribution (figure [3]) is due to the large scatter in 
the individual measured values, reflecting the scatter in 
the V — L relation. These figures visually demonstrate 
the need to average many magnification measurements 
in order to obtain a result precise enough to compare to 
theory/models. 

The transverse distance between each quasar and the 
MaxBCG cluster closest on the plane of the sky is shown 
in figure [5] panel (a) shows the distance in the logarithm 
of kiloparsecs at the redshift of the closest cluster, while 
panel (b) shows the distance as a fraction of i?2oo of the 
closest cluster. 

5.5. Magnification profile of the MaxBCG Clusters 

For each quasar, we calculate the magnification we ex- 
pect due to clusters in the MaxBCG catalog, as described 
in section [5^ We then group the quasars in logarithmi- 
cally spaced bins given by the quasar's distance from the 
closest cluster, in units of i?200 of that cluster. The errors 
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Figure 4. Logarithm of the tensing magnification of quasars as- 
suming NFW profiles for iVIaxBCG galaxy clusters. The thick-lined 
histogram includes one entry for each of the 3573 quasars which 
overlap the iVIaxBCG area. The thin-lined histogram shows the 
same calculated magnifications, but includes one entry for each of 
the 57,359 quasar measurement pairs and is therefore comparable 
to figure [31 

on the measurement are determined through bootstrap 
resampling. In particular, for each bin 1000 mean mag- 
nifications are calculated using random sampling with 
replacement such that the randomized set includes the 
same number of measurements as the original bin. The 
errors are taken to be the la (15.9 and 84.1 percentile) 
values of the 1000 mean magnifications. Because the er- 
rors on the magnifications are in part caused by system- 
atic uncertainties that are not currently understood, we 
prefer to derive the errors empirically in this way using 
the properties of the distribution itself. 

Figure [Hi shows the mean measured /i' (points with er- 
ror bars) and the expected magnification (points without 
errors, connected by lines), assuming NFW haloes for the 
MaxBCG catalog's clusters as described above. While 
the magnification is significant and does fall with radius, 
we measure a systematically different magnification pro- 
file from that expected. The measured values compare 
to the expected magnifications with reduced = 1-84. 

The MaxBCG catalog's cluster positions are given as 
the location of the brightest cluster galaxy (BCG). This 
position does not always coincide with the center of the 
dark matter halo. Errors in the cluster central coordi- 
nates cause the average radial profile of the MaxBCG 
clusters (as me asured from the repor ted centers) to dif- 
fer from NFW. [Johnston et"all ([2007D studied this effect 
using weak lensing shear measurements and N-body sim- 
ulations, and found a Gaussian distribution of positional 
offsets between the true and measured centers, with a 
richness-dependent probability that the cluster is cor- 
rectly centered. To determine the effects of such offsets 
on our exp ected magnificat i ons w e implement the pre- 
scription of [Johnston et ahl ([2007[ ). offsetting the center 
of a subset of the clusters in a random manner. The 
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Figure 5. Distance, transverse to the line of sight at the redshift of 
the cluster, between the quasars and the clusters which are closest 
on the plane of the sky. (a) in units of log kiloparsecs at the redshift 
of the cluster, (b) in units of iJ200 of the closest cluster. 
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Figure 6. Measured /i' (points with error bars) and expected fi 
(points connected by lines), versus distance from the nearest cluster 
in units of R200, assuming NFW profiles truncated at 3 X i?200- 
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results are shown in figure [71 and provide an improved 
fit to the data, with reduced — 1-17. The agree- 
ment between the data and predictions may be further 
improved by modelling the effects of large scale struc- 
ture outside the redshift range of the MaxBCG catalog; 
however, these effects are subdominant and beyond the 
scope of this paper. 

Using the errors on the binned measurements deter- 
mined through bootstrap resampling, we can estimate 
the error on each magnification measurement. Figure 
[5] shows the logarithm of the error on each bin versus 
the number of measurements in that bin, taken from the 
data in figure [71 with a best-fit line superimposed. The 
trend falls as Vn, and at the level of our most populated 



Figure 7. Measured /i' (points with error bars) and expected fi 
(points connected by lines), versus distance from the nearest cluster 
in units of R20O1 assuming NFW pr ofiles truncated at 3 X i?200 and 
with centers off'set as prescribed in ljohnston et al] ||2007I 1. 

bin (with 21,482 measurement pairs) shows no sign of a 
systematic floor. The best-fit line extrapolates to imply 
that each log(/j,') measurement, made using one magni- 
tude measurement pair, has an error of roughly 0.4. For a 
magnification of 1 . 1 , this corresponds roughly to a signal 
to noise of 1 for each magnification measurement. 

6. DISCUSSION AND CONCLUSIONS 

We introduce a new technique to measure the effects 
of gravitational lensing by large scale structure. Type I 
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Figure 8. Logarithm of thie error on eacti bin versus the number 
of measurements in that bin, taken from the data in figure O with 
a best-fit line superimposed. The trend falls as \/7V, implying an 
error on log(/i') of 0.04 for a magnification measurement using one 
pair of magnitude measurements. 

quasars on average follow a well-constrained variability- 
luminosity relation; lensing magnification causes a sig- 
nificant shift in this relation. We describe the magnifi- 
cation measurement procedure, and use it to determine 
the magnification of background quasars due to lensing 
by the galaxy clusters in the MaxBCG catalog. 

The measured magnification is significant, and consis- 
tent with the signal expected from the catalog. In partic- 
ular, the agreement between the data and expectations is 
improved when the miscent ering of the clusters is taken 
into account, as described in lJohnston et al.l()2007D . com- 
pared to when the clusters are assumed to be NFW pro- 
files centered exactly on their catalog positions. The lens- 
ing effects of the MaxBCG catalog have been measured 
very precisely using weak lensing shear, for example in 
iSheldon et all (|2009[ ). Ahhough the magnification mea- 
surements presented here are not competitive with such 
results, they serve to show that quasar variability can 
be used to study lensing magnification, and are the first 
measurements using this new technique. 

Using quasar variability to measure gravitational mag- 
nification is a very different approach from the common 
method of using galaxy shapes to measure gravitational 
shear. Precision shape measurements of distant galax- 
ies require very good quality data with faint flux limits 
and small point spread functions. Furthermore, the com- 
plex intrinsic shapes of galaxies make the shear measure- 
ment nontrivial. Measuring the variability amplitudes of 
quasars is much more simple, as these objects are often 
bright and they typically vary by amounts larger than the 
photometric errors of surveys such as Palomar-QUEST 
or SDSS. We see the error on the measured magnification 
to be roughly given by 

log(aiog^) = -0.51og(7VA^)-0.4 (10) 



where iVAm is the number of Am measurements. This 
can be compared to the error on weak lensing mea- 
surements of the reduced shear g, which scale as ag = 

CTe X (1 — |5p)/\/-^gai ~ cre/-y/A'"gai for moderate shears. 
The intrinsic scatter i n the galaxy shapes is ~ 0.35 
(jSchneider et al.|[2004[ ). and iVgai is the number of mea- 
sured galaxies. For magnification by a galaxy cluster 
with an NFW profile parameterized by M200 = 10^* and 
concentration c = 6, a region with moderate magnifi- 
cation fi = 1.05 would have reduced shear g = 0.035. 
In this regime, 100 galaxies and 100 Am measurements 
would yield signal to noise S/Nfj_ — 10 and S/Ng — 
1. A region of higher magnification /i = 1.25 and g = 
0.083 would yield again S'/iV^ = 10, with S/Ng ^ 2.4. 
This lensing strength is at the upper end of the range 
for which weak lensing techniques have b een rigorously 
teste d, for example by the STEP2 project (jMassev et al.l 
l2007f ). Stronger distortions are not straightforward to 
measure, particula rly using the common KSB method 
aiser et al.lll995f) . The variability-based quasar mag- 
nification analysis, on the other hand, is equally applica- 
ble to all lensing regimes. 

The fact that the magnification measurements each 
have higher signal to noise than the shear measurements 
is offset by the relative sparsity of quasars compared 
to inactive galaxies. In the current SDSS spectroscopic 
sample, which reaches depths of r ~ 19.5, there are 
about 15 quasars per square degree. In data down to 
r ^ 20.5 we expect ro ughly 40 quasars per square degree 
(|Richards et al.|[200^ . This depth will soon be available 
over large areas of sky. For example, BOSS is expected to 
double the number of spectroscopic quasars in the SDSS 
Survey over the next several years, over a footprint of 
10,000 square degree^ Shear analyses are usually per- 
formed over much smaller areas of high quality data. For 
example, results from the CFHT Wide Survey use 50,000 
galaxies per square degree, down to a magnitu de of i ~ 
24.5, with a total area of 34 square degrees (|Fu et al.l 
I2OO8I ). In order to reach the same signal to noise for 
moderate lensing strengths oi n — 1.05 and g — 0.035, 
assuming 40 quasars and 50,000 galaxies per square de- 
gree, we would need fewer than 4 magnitude measure- 
ments on average for each quasar (assuming that the 
measurements give uncorrelated information about the 
magnification; more work is required to investigate this 
systematic complication) . Assuming shallower data with 
15 quasars per square degree, we would require 6 mag- 
nitude measurements of each quasar, on average, over 
just 34 square degrees. Alternatively, the same signal to 
noise could be achieved with just 2 measurements per 
quasar over an area of 1100 square degrees. These lev- 
els of data are easily obtainable with surveys such as 
BOSS and Pan-STARRS, which will provide >10,000 
square degrees with sufficient data quality for magnifi- 
cation measurements. This is in sharp contrast with the 
70 square degrees which will be available from Pan- 
STARRS Medium Deep fields with the quality assumed 
here for shear measurements. 

This work introduces the technique of using variabil- 
ity to measure the lensing magnification of quasars, and 
measures the magnification of 3573 quasars, with magni- 

^ http://www.sdss3.org/ 



12 



tudes down to SDSS r ^ 19.5, due to Icnsing by known 
galaxy clusters along the line of sight. Fiitiire large-area 
sky surveys will be able to apply this new technique to 
measure lensing magnification over larger volumes, and 
improve our understanding of large scale structure. 
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